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SECTION I 


1.0 SUMMARY 

A method haa been developed to estimate the weight and major 
envelope dimensions of small aircraft propulsion gas turbine 
engines. The computerized method, called WATE-S (Weight Analysis 
of Turbine Engines - Small) is a derivative of the Boeing WATE-2 
computer code (ref. 1), As in WATE-2, WATE-S determines the 
weight of each major component in the engine including com- 
pressors, burners, turbines, heat exchangers, nozzles, pro- 
pellers, and accessories. A preliminary design approach is used 
where the stress levels, maximum pressures and temperatures, 
material properties, geometry, stage loading, hub/tip radius 
ratio, and mechanical overspeed are used to determine the com- 
ponent weights and dimensions. 

A relatively high level of detail was found necessary in 
order to obtain a total engine weight within the required 
±10-percent accuracy. Component weight data from Garrett engines 
was used as a data base to develop the methods required for small 
gas turbine engine components. The small engine weight method is 
probably valid for other manufacturers, however, this has not been 
verified. The accuracy of the method is generally better than 
±10 percent, on the order of ±5 percent. This accuracy was 
verified by applying the method to four Garrett propulsion 
engines, some of which were in the original data base. Engines 
used in the validation study were jointly selected by NASA and 
Garrett . 
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SECTION II 


2.0 INTRODUCTION 

Aircraft propulsion system stuaies are frequently conducted 
by industry and government. These studies may encompass a wide 
variety of engine concepts ranging from relatively simple turbo- 
prop and turbofan engines to complicated variable cycle engines. 
The aerospace industry in general has acquired an adequate com- 
puter capability to evaluate the thermodynamic performance of 
these diverse engine concepts; however, until recently, accurate 
methods of estimating engine weight and dimensions were not read- 
ily available. 

The Boeing WATE-2 (Weight Analysis of Turbine Engines) com- 
puter program was written primarily to estimate the weight and 
major envelope dimensions of large, axial-flow turbofan and turbo- 
jet propulsion engines. A preliminary design approach is used 
where the stress levels, maximum temperatures and pressures, mate- 
rial properties, geometry, stage loading, hub/tlp rad? us ratio, 
and shaft mechanical overspeed are used to determine component 
weights and dimensions. 

Overall, the basic approach of the Boeing WATE-2 program is 
valid for small gas turbine engines; however, major modifications 
to some of the component modules were required to Improve the 
accuracy and flexibility of this program for small gas turbine 
engines . 
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A computer program apeoifically written to patimate the 
weight and dlmenalona of amall gas turbine engines was derived 
from the Garrett version of the Boeing WATB-2 computer program. 
This program will be referred to as WATE-S (Weight Analysis of 
Turbin<.} Engines - Small) . It was configured on a stand-alone 
basis, the required input being obtainable from most general pur- 
pose thermodynamic engine cycle performance computer programs. A 
small gas turbine engine as defined for WATE-S suitability is out- 
lined below: 

o Conventional and unconventional Brayton Cycles (regen- 

eration, intercooling, reheat, afterburning, combina- 
tions) 

o Separate or mixed-flow turbofans 

o Core-corrected airflows from 0.23 to 4.54 Kg/s (0.5 to 
10.0 Ibra/s) 

o Compressor pressure ratios from 5 to 25 

o Turbine rotor inlet temperature from 1089“K (1960“R) to 

1922»K (3460®R) 

o Propellers and propfans suitable for aircraft flight 

speeds up to Mach 0,8 

o Axial and radial flow tucbomachinery components with 

input specified staging and shafting arrangements 

o Geared or ungeared propulsor drive. 

WATE-S determines engine weight and dimensions by summing the 
weights and dimensions of each major component including pro- 
pellers, fans, compressors, combustors, turbines, and nozzles. 


Each Individual component specified in the input is Individually 
modeled using fundamental physical relationships. The specific 
tasks accomplished are summarized in the following sections, 

2,1 Component ModuJis Requiring Modifications 

Fans and Axial Conipressors - The WATE-2 axial compressor model 
assumes the general design pressure-ratio corrected tip-speed 
relationship is applicable to fans as well as low- and high- 
pressure compressors. Garrett experience has demonstrated that 
this assumption is not valid for most small and intermediate tur- 
bofan engines. Separate fan accounting is provided in WATE-S as 
well as a design pressure-ratio corrected tip-speed relationship 
suitable for small axial compressors. 

* 

Centrifugal Compressors and Radial Turbines - The WATE-2 centrif- 
ugal compressor and radial turbine models are not consistent with 
the desired component preliminary design philosophy. These mod- 
ules have been completely rewritten so that radial flow component 
weights and dimensions ate determined on the basis of fundamental 
physical relationships similar to those used for axial flow compo- 
nents . 

Reverse-Flow Co.. . justora - The WATE-2 program’s combustor logic is 
not applicable to reverse-flow combustors. The mean burner diam- 
eter is used as the geometric design parameter. For reverse-flow 
burners, where the burner is positioned around the HP turbine, the 
inner diameter has more geometric significance. The ability to 
use either of these diameters as the geometric design parameter Is 
incorporated in WATE-S. 

Engine Accessories - Engine accessory weights are currently lumped 
together as a percentage of bare engine weight. Small gas turbine 
engine accessory weight data including the starter/generator , 
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battery, tuel pumps, oil pumpe, propeller pitch control, fuel con- 
trol, etc., was reviewed to determine if a more detailed weight 
model waa warranted. A new accessory weight model has been devel- 
oped and Included in WATE-S. 

2.2 New Component Modules 

Gearboxes - The WATE-2 program gearbox weight accounting is spread 
through a number of subroutines, A different model is used for 
axial and centrifugal compressor gearboxes, and no method is 
available for linking a gearbox to a turbine for geared fan or 
propeller f tee-turbine engines. A new gearbox component module 
was necessary to provide the required flexibility for small gas 
turbine engines. In this way, gearboxes may be connected to any 
rotating component desired and consistent gearbox weights esti- 
mated , 

Propeller and Propfans - Methods of estimating propeller and prop- 
fan weight and dimensions are not included in WATE-2. A new pro- 
peller component module based on Hamilton Standard propeller data 
has been developed and incorporated in WATE-S, 

2.3 General Modifications 

In addition to the above, the printer plot option of WATE-2 
has been made operational for small gas turbine engines. Centrif- 
ugal compressors, reverse-flow combustors, and radial turbine ele- 
ments were added to the axial flow component plot capabilities. A 
general review of the many geometric assumptions used in WATE-2 
for large turbomachinery components was completed to determine 
their suitability to small gas turbine engine components. 


A default engine configuration is incorporated in WATE-S. 
The default engine is a two spool turbofan engine consisting of a 


fan, axl-Qentrlfugal hlgh"preaaur« 

PceBsura turbine, and axial fan turhi «‘«l h» 9 h- 

hava been uaed for oomputet proaram oh"*'t f»f*'«nce englnea 

accuracy of WATB-s In el, 1 '">U«atlon. The 

Oimenalona la «o percent or bLtT''"-^i^T 

by analyals of the reference enplnea’. '’®" Oemonatrated 




SECTION III 


3.0 METHOD OP ANALVSIS 


The methodB of analysis used in WATE-S are based on a preli- 
minary design approach where the aerodynamic and mechanical design 
variables are taken into account, in the compressor, for example, 
rotor blade weight is calculated as a function of the specified 
geometric parameters. Blade centrifugal stress is then found, and 
the disk weight that will support the blade rotational force is 
determined. This type of preliminary design approach was also 
used for the other components. 

The WATE-S method is intended to estimate the weight of a 
given engine design; it will not design an engine. This function 
must be performed external to the program. WATE-S utilizes com- 
ponent state conditions which are generated in an external engine 
thermodynamic cycle computer program. 

In the normal use of WATE-S, the desired engine cycle is sim- 
ulated at the engine design point. The user of WATE-S must be 
cognizant of other conditions In the flight envelope where maximum 
component temperature, work, speed, or flow occur. If these con- 
ditions are greater than the design values, they can size the com- 
ponent and have a significant impact on the component weight. 
WATE-S allows input of scalars to account for these off-design 
conditions. 


A more accurate weight estimate can be achieved by developing 
an array of engine cycle data over the engine operating envelope. 
The WATE-S program will scan the input engine cycle data and 
select the maximum conditions for each component. 


The engine oyalo nlmulatlon may mao coroponenta that are 
required mathematlaaUy , but are not dealred In the engine 
weight. Those can be selectively eliminated by the use oC the 
component weight scalars. 

3 ^ Compon ent weight and Dimen sion Metho d olog y 

The methods o£ analysis described for each component In the 
following sections have been developed to achieve an overall 
engine accuracy of ±10 percent. Since the rotating components 
comprise the major part of the total engine weight, considerable 
detail was necessary in order to achieve the accuracy goal. Not" 
mal program users may not have sufficient knowledge to adequately 
define all of the necessary inputs; however, typical values are 
given in the User's Guide, Section 4. 

3.1.1 Axi al Fans a n d Compressors 

The procedure used for axial fan and compressor weight pre- 
diction Is a stage-by-stage aeromechanical design as illustrated 
in Figure 1. Rotor blade volume and weight are determined; then, 
blade centrifugal stress, disk stress, and disk weight are calcu- 
lated. Connecting hardware, stator blades, and cases are then 
estimated and summed to give the total component weight. The fol- 
lowing input data is necessary: 

o The number of compressor stages. Stage work is held 
constant for multistage axial compressors. Alterna- 
tively, a maximum first-stage pressure ratio, which 
reflects the design approach and technology level, may 
be input 

o Inlet and exit Mach numbers of the compressor 
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THERMODYNAMIC MATERIAL STRUCIU 

GAS PATH CHARACTERISTICS REQUIRE 

CHARACTERISTICS 



Compressor Weight Method. 







o 


Inlet hub/tlp 


cadlua ratio of: the Clrat-atage blade 


o CoirptMSOC aealsn moao, l.e., oonatant maanUne, oon- 
Btant hub, or oonatant tip diameter 

o ECteotlve aenslty o£ blade matetlali daSlned as the 
ratio of total blade weight and volume 

o Maximum Inlet and exit temperatucea , If not at design 

o Aspect ratios for the first and the last stage blades 

o RPWmax/*^^“des overapeed factor 

o Blade tip solidity, ratio of blade tip chord and spacing 
o Density of disk material 

o Blade taper ratio 

o Blade volume factott ratio of blade volume and total 
annulus volume. 

The total enthalpy change foe the component is available from 
the engine cycle data. Equal work for each stage Is assumed, anu 
the number of stages required Is found by Iteration until the 
first-stage preseure ratio Is equal to or less than the speoltie 
maximum. When the number of stages Is specified, first-stage 
pressure ratio la calculated based on the equal work per stage 
assumption and the maximum allowable pressure ratio Input Is 
Ignored. Shaft speed Is determined from the axial fan and com- 
presaor tip speed oortelatlone given In Figure 2. These com- 
pressor tip speed correlations are baaed on current (1982) tech- 
nology levels. The rpm of additional compressors driven on the 
same ahaft will be set by the first upstream compressor. 
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In the event that an existing compressor la being weighed, or 
the shaft speed Is already known, a speed scalar can be applied to 
the shaft speed found from Figure 2 to achieve the desired value. 
The speed scalar can also be used to modify the estimate of shaft 
speed for high-pressure compressors, where the Inlet temperature 
significantly affects the pressure ratio capability, or for 
external Iterations of an engine design. 


The first-stage flow area Is determined by the specified Mach 
number and the corrected Inlet airflow from the cycle data. Inner 
and outer diameters of the flow path are calculated from the spec- 
ified radius ratio: 


D = 

^ V7rll-(h/t )^] 


( 1 ) 


D 


h 



( 2 ) 


Ccanpressor RPM is determined by dividing the tip speed (found from 
Figure 2) by the product of tt and the tip diameter (D^) . 


Stage length is found In the following manner. Axial blade 
chord (C), as shown In Figure 3, Is the quotient of the blade 
height and aspect ratio: 
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The stator length was assumed to be equal to the rotor length 
(or blade chord), and 17 percent of the rotor length required for 
clearance between rotor and stator and the same clearance between 
the stator and the next rotor. 


The total number of blades is calculated from the specified 
tip solidity (C/S) and aspect ratio (AR) , and previously calc-'i~ 
lated tip diameter (D^) and blade height (hg) : 


N 


B 


7T 


c/< 


AR 


( 4 ) 


This value is truncated to an integer number of blades, and the 
same value is used for the stator. 


The total volume of metal in the compressor blades is then 
calculated from; 


V, 


‘B 


B 


AR 


?■ 


(5) 


where K is a volume factor which accounts for flrtree mount vol- 
ume, taper ratio, and thickness-to-chord variations in typical 
blades. For the data base engines, K was found to be 0.06 for fan 
blades and 0.12 for compressor blades with hub/tlp radius ratios 
less than 0.8. For compressor blades with hub/tip radius ratios 
greater than 0.8; 


K » 0.12 + 0.04 (h/t - 0.8) 


( 6 ) 


The rotating blade weight of each stage is determined from 
the blade volume and material density. Material density automati- 
cally changes from titanium to steel when the stage Inlet tempera- 
ture exceeds a specified maximum, normally 644°K (1160**R). Stage 
temperature Is determined from the engine cycle data, the equal 
work assumption, and the specified over-temperature ratio. Stator 
blade weight and dimensions are assumed equal to the rotor blades 
and Include the inner shroud. 



Tlw maximum blarto centrifugal atrean occura at the hi arte root 
„nrt in a function o£ tip apoed, blade height, taper ratio, and 
material denalty. Expreaaert In terma o£ the nondlmenalonal blade 
geometry, the equation for blade toot oentrltugal atroaa la. 


lOOpU^' 




TR-l 

’TT' 


(1 - h/t) (1 + 3 h/t) 


where U. is the product of the design tip speed and the overspeed 
factor, 

material Is assumed to be titanium or steel forgings with ultimate 
strengths of 82.7 KN/cm^ and 110.3 KS/cm^ respectively. The rim 
of the disk is assumed to be 10 percent of the blade hub radius or 
1.91 cm, whichever Is greater. The remaining disk volume Is 
assumed to be a trapetoldal section rotated about the axis. The 
thickness of the trapezoid at Its outer radius is sized fo 
75 percent of the ultimate strength. Thickness of the trapezoid 
at the Inner radius is based on 50 percent of ultimate strengt 
for the tangential stress. This value was selected based on an 
experimentally determined burst speed margin. 

If the average tangential stress can be satisfied with a con- 
Slant thickness disk, the bore radius Is Increased until the 
lowest weight disk is achieved that just meets the design crl- 

ter ia. 

Disk stresses are calculated based on the total blade force 
acting on the outer circumference of the disk, and are therefore 
sensitive to changes in blade aspect ratio and solidity, as well 
as blade centrifugal stress and hub radius. 

Dimensions of succeeding stages are based on the design mode 
nuloctPd (constant moan, tip, or hub radius, In the constant 
moan lino method, for oxamplo, the moan radius Is ^ 

moan flow area of the first stags. This moan radius is hold 


constant for subsequent stages. Corrected airflow at the entry of 
each stage is determined from the calculated state conditions 
derived from the equal work per stage assumption. Stage Inlet 
Mach number Is assumed to vary proportionally to the number of 
stages when different Inlet and exit Mach numbers are specified. 
Tip speeds for the downstream stages ace then calculated from the 
stage dimensions and shaft speed. Blade aspect ratio Is deter> 
mined by assuming a proportional change for each stage If Inlet 
and exit aspect ratios are different. 

Figure 4 illustrates the stage coupling method that was dev- 
eloped by Boeing and used In WATE-S. The spacer, nuts, and bolts 
ace assumed to be steel, with the spacer being a 0,191-cm thick 
cylinder located at 75 percent of blade hub radius. The connec- 
ting hardware weight is estimated by the following equation: 

Whw = 27 t (0.75 . R^) . 0.191 , .p (8) 

where R^^ is the blade hub radius, is the stage length, and p 
is the material density. 

The outer case is the last item of weight included in the 
compressor weight buildup. Average case thickness in the data 
base engines was 0.254-cm equivalent thickness, including 
fasteners and flanges. Case weight is calculated stage by stage, 
and the same material used in the disk is also assumed for the 
case : 


Wc = . Lg,pQ . 0.254 .p (9) 

where 1s the stage tip diameter. 

Total stage weight Is the sum of the rotor blades, stators, 
disk, connecting hardware, and case. Stage weights are summed to 
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give the total component weight. The sum of all stage lengths 
gives the total compressoi: length, inlet guide vanes are not 
Included in the compressor weight , hut can be accounted for with a 
frame <see Section 3.1*9). 

Rotational Inertia of the compressor la determined by finding 
the inertia of each stage of the component. It la assumed that 
blades have uniform weight/area and the disk is of uniform thick- 
ness. This method had been applied to several existing engines 
and the results showed good agreement. The following equations 
are used for the blades and disk Inertias: 




[r(^ “h*)] "b 


( 10 ) 


T « 2 

" 2g * ^h 


( 11 ) 


The total stage Inertia is the sum of blade and disk Inertias 


I ® I + I 
^Stg D 


( 12 ) 


The component Inertia is the sum of the Inertia of each stage. 


3.1.2 Rotating Splitter 

A rotating splitter, see Figure 5, Is a circumferential 
separator of two flows within the same component. These flows 
normally have different pressures and temperatures, and the split- 
ter must perform a sealing function. Stages that Incorporate 
rotating splitters are treated the same as compressors; a rotor- 
stator pair comprises one stage, stator weight and size are 
assumed to be the same as the rotor blade, and rotor-stator spac- 
ing is 17 percent of rotor length. 


13 


OF POOR QUALITY 



19 



The rotating spllttec adds wel9ht to the blade and Ineroaaee 
the centrifugal blade force. Consequently, the disk must be 
heavier to carry the added load. Splitter weight per blade is 
estimated byi 


”SPL 



0.10 



(13) 


where C 1s the blade chord found from Equation (3) and Rgpj^ la the 
radial location of the splitter. Thickness of the splitter Is 
assumed to be 10 percent of the chord, however, this choice was 
based on only one engine, the General Electric CJ805-23. The 
CJ805-23 aft-fan blade has a rotating splitter which has a box 
section. The solid equivalent thickness of L..e hollow box was 
approximately 10 percent of the chord. A more accurate estimate 
could be made by actual design of the cantilevered platform to the 
desired deflection ana/or stress levels. 


The centrifugal force contribution of the rotating splitter 
is : 


SPL 


Mrw" 


W 


—SPL 
lO'Og 


2v. 


'SPL 


RPM 

60. Q-' 


(14) 


where RPM is the shaft speed determined in the same manner as for a 
compressor . 


Blade root centrifugal stress (Equation 7) is increased by 
the amount 


!§PL 




:SPL 

. t/c 


(15) 


where t/c is the thickness/chord ratio of the blade (which is 
assumed to be 10 percent). Disk weight is determined with the 
increased stress level using the procedures described for axial 
fans and compressors (Section 3.1.1). 
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Shaft speed determination (as described In Section 3.1.1) 
assumes that blado'-root stress 1s subcrltlcal. Use of a rotating 
splitter will cause the blade centrifugal stress to Increase slg- 
niflcantly, and the WATE-S output should be inspected to determine 
whether or not the stress level Is acceptable. Reduction of shaft 
speed may be required to reduce the stress level. 

If shaft speed is decreased, additional compressor stages may 
be required to accomplish the same work. Alternatively, the 
radius ratio of the compressor can be Increased to restore work 
capacity (due to higher tip speed). Disk weight of each stage 
will increase for this compromise, however. The final choice must 
be Iterated external to WATE-S and may depend on whether or not 
the flow path Is reasonably well matched to connecting components. 

These secondary effects may have a much larger impact on the 
engine weight than the weight of the splitter material, and they 
should not be ignored. 

3.1.3 Axial Turbines 

The methods used for axial turbines are generally the same as 
those previously described for axial compressors. The following 
input data is required: 

o The number of stages or the maximum mean diameter of the 
first stage 

o Inlet Mach number (axial) of the first stage, and exit 
Mach number (axial) of the last stage 

o Rotor blade aspect ratios of the first and last stages 

o Blade tip solidity 
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o Rofarenae tUak atreaa, 0,2 percent yield point of the 
materiel eoleoted 

o Cooling Indicator, to modify the blade volume oalcula- 
tlon for cooling passages 

o Turbine design mode) l,e., constant hub diameter, con-^ 
stent mean line, or constant tip diameter 

o Shaft overspeed factor, 

o Mean work coefficient, ■ 9--~- 

m 

where 4H/N » enthalpy change pet stage 
and *■ mean blade speed 

o Blade material density 

o Blade taper ratio 

o Blade volume factor; ratio of blade volume and total 
annulus volume. 

Two alternative procedures may be used to size the turbine: 
(1) specify the number of stages, and assuming equal work per 
stage, calculate the mean diameter of the first stage based on the 
Input turbine mean work coefficient; or (2) specify the maximum 
mean diameter of the first stage, and assuming equal work per 
stage, calculate the number of stages required to satisfy the 
input turbine mean work coefficient. Total component work and 
state conditions are taken from the engine cycle data. 

The flow area required to pass the corrected airflow at the 
specified Mach number is calculated at the inlet of each stage. 
For the first stage, the hub and tip radii are calculated from 
this area and the previously determined turbine mean radius. Sub- 
sequent stage dimensions are calculated baaed on the design mode 


selected (ponetant mean, hub^ oi: tip diameter). When the tlret- 
atayp Inlet Mach number and aspect ratio are different from the 
last stage values ^ a proportional change is assumed for the Inter- 
mediate stages. 

Kotor blade chord and number of blades are determined by the 
same methods used for axial compressors. Blade volume Is also 
determined by the same method^ Equation (5), except that K ■ 0.155 
for high-pressure turbines and K ■> 0.195 for low-pressure tur- 
bines. When the blade is a cooled blade, the calculated volume is 
reduced 20 percent to compensate for coollng-alr passages. Blade 
weight is then found from the specified material density. Blade 
root centrifugal stress is calculated by Equation (7). The disks 
are calculated by the preliminary design procedure already discus- 
sed for axial compressors. 

Each stage of the turbine is treated as a stator-rotor pair 
(as opposed to rotor-stator pair in the compressor) . Stator 
blades are assumed to have the same number and volume of material 
as the rotor blades. The stator weight is calculated by Equation 
(5), with K B 0.155 for high-pressure turbines and K » 0.195 for 
low-pressure turbines. Stator-rotor spacing is the same as com- 
pressors, 17 percent of the rotor length. Stator AR is taken to 
bo 83 percent of the rotor blade AR. 

Connecting hardware and case weight are also determined by 
the same methods used in the compressors. The total weight and 
length of the turbine component is the sum of disk, blade, stator, 
connecting hardware, and case. No exit guide vanes (EGV) are 
assumed in the turbine component. EGV's, if required, can be con- 
sidered a part of the exit frame weight. Rotational inertia is 
determined In the same manner as the compressor. 
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Th« procedure used Cor cantrlCugeX compreesore Is a prelim- 
inary design approach similar in principle to that used Cor axial 
compressors. Blade volumes and weights are determined and a disk 
is sized to support the blades Cor allowable stress levels based 
on inlet and exit Clow path temperatures. The Collowing inputs 
are required Cor centrifugal compressorat 

o inlet Mach number 

o Maximum first-stage pressure ratio, which reflects the 
design approach and technology level. Stage work for 
this stage will be held constant for the second stage. 
Alternatively, the number of stages can be specified. 

o Inlet hub/tip radius ratio 

o Number of blades and splitters 

o Exit Mach number 

o ^^^roax/'^^^des over speed factor 

inlet/^^Pexit 

Only a one- or two-stage centrifugal compressor may be used 
In a WATE-S engine simulation. The total enthalpy change for the 
component is available from the engine cycle data. Equal work for 
each stage is assumed, and the number of stages is calculated to 
satisfy the maximum first-stage pressure ratio requirement. When 
the number of stages is input, first-stage pressure ratio la cal- 
culated based on equal work per stage and the input maximum pres- 
sure ratio is Ignored. 



Tho Inlot fli)w art'rt i« 4 otormino <1 fiom tlu’ Hpoctflert Mnch 
numluM anti Inlat. ourroctejl flow from th«> onglne cycle dntn. The 
inlot hub |l*j||> «n,’l uhroud rtt< 1 ll are calculated from the 

apoclfliHi Inlet/hub tip radlue ratio (see Figure 6 ) . 


o»- 


Kigur*' . Tontr i fugal ComprcHSor/Hadial Turbine Schematic. 


. 

0 U-(hA)^l 


"ill ’ "/t . K,js (17) 

The exit tip radiun (K^) la then calculated from the input 

*’^’*'inlot'^*^ **Vxlt ratio. The exit flow area ia detorralnod 

from the apocified exit Mach number and corrected flow from the 
engine cycle data. The compreasor exit blade length (in.) la cal- 
culated from this area and R. i 

t 

\xlt 


OF pOt’K 



The ratio of axial length to radial height of the blade at the 
shroud t assumed to be 2., and the flow path is assumed to 
be comprised of two ellipses. Blade volume is computed assuming 
the blade thickness at the shroud is 0.076 cm and the thickness at 
the hub results in 15-percent blockage; 


^hub 


0.15 (27rRiH^ 


(19) 


where is the number of blades. The compressor material and 
allowable blade and disk stress levels are selected based on the 
maximum temperature conditions for the compressor. 


The compressor speed is determined from the current (1982) 
technology design pressure-ratio corrected tip-speed correlation 
illustrated in Figure 7. Once the material is selected and the 
rpm is calculated, the blade weight, stress, and natural fre- 
quencies are calculated. 



FIRST-STAGE CENTRIFUGAL COMPRESSOR PRESSURE RATIO 

Figure 7. Centrifugal Compressor Tip-Speed Correlation. 
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The GompreHSor disk la sized by an Iterative preliminary 
design procedure to support the blades and minimize the stage 
weight, while maintaining the maximum allowable stress level at 
the bore. This Is done by varying the bore radius and backside 
geometry o£ the disk. 

Total stage weight is the sum of the blade, disk, shroud, and 
diffuser weights; and crossover duct weight if a two-stage centrl'- 
fugal compressor is modeled (see Section 3.1.6, Equation 32). 

R + R ^ 

“shroud - 

“dlf«UBet “ 0'2845p (21) 

The diffuser radius is assumed to be 1.6 times the exit tip radius 
(R.) and the length equal to 7 percent of its diameter. 

The total stage inertia is estimated from the following 
expression: 


Icent " 0-1086pR^ <22) 

where (> is the material density selected. 

3.1.5 Radial T u rbines 

The method described for centrifugal compressors is generally 
the same for radial turbines. The input data required is; 

o Inlet Mach number 

9AH/n 

o Turbine work coefficient, ■ 2 " 
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o 


Number of blades and splitters 


o Exit Mach number 

o overspeed factor. 

Only a single-stage radial turbine Is allowed. The total enthalpy 
change for the turbine is known from the engine cycle data. The 
Inlet tip radius (R^) of the turbine is calculated from the input 
turbine work coefficient and the previously calculated shaft 
speed. The inlet flow area is determined from the specified Mach 
number and the corrected flow from the engine cycle data. The 
inlet blade length is calculated from: 

BL = ~~ 

27rRt ( 23 ) 

The exit shroud radius turbine is assumed to be 

70 percent of the inlet radius (R^) and the exit hub radius 
is calculated knowing the exit corrected flow and Mach number. 

The ratio of axial length to radial height of the blade is assumed 
to be 2.1, and the flow path is assumed to be comprised of two 
ellipses. Blade volume is computed the same as for centrifugal 
compressors. The turbine material is assumed to be a nickel-based 
superalloy with a density of 0.00792 kg/cm^ and allowable stress 
levels for the blade and disk of 51.7 KN/cra^ and 86.2 KN/cm^, 
respectively. Disk weight is calculated the same as for the cen- 
trifugal compressors. 

Total component weight is the sum of the blade, disk, and 
shroud weights. 
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The total inertia of the 
following expression: 
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0.1755^ 



(25) 


radial turbine is estimated from the 


^rad turb 


0.1687/>R^5 


(26) 


3,1.6 Ducts 


The design criteria used to size the ducts in WATE-S is 
assumed to be the Internal pressure. Structural loads cannot be 
addressed in the preliminary design process where WATE-S is 
intended to be used. Figure 8 illustrates the typical duct geome- 
try assumed. Generally » the outer surface of the OD duct wall is 
exposed to ambient pressure or fan duct pressure. The inner sur- 
face of the ID duct wall can be subjected to fan pressure, LP 
compressor pressure, or HP compressor exit pressure, etc. For 
these reasons, ambient pressure is assumed for the calculations of 
duct wall thickness, this assumption results in conservative 
weight estimates. 



Figure 8. Duct Schematic. 
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The equation tor atresa on a longitudinal aectlon of a thin- 
walled cylinder aubjected to internal preaaure is (ref. 2): 


ir = 


PD 

2t 


(27) 


or solving for the minimum thickness, 


t 


_ PD 
min ' 2^r 


(28) 


where rr is the allowable stress level, P the internal total pres- 
sure, and D the duct wall diameter. A minimum gage thickness of 
0.127 cm is assumed for all ducts. The assumed materials are 
titanium {(J'= 34.5 KN/cro^) below 644*K (1160“R) and stainless 
steel (cr» 48.3 KN/cm^) above 644“K (1160®R). The material is 
selected based on the total temperature of the duct airflow. The 
weight is calculated as a function of the duct length (L) , the 
inner diameter (D^), the outer diameter (D^) , and the calculated 
wall Ihloknessea and 


«duet '°o‘ 


min ,o 


^i^min,i^ 


(29) 


Duct Mach number is specified as an input, and corrected air- 
flow is determined from the engine cycle data. The inner and 
outer diameter are determined as a function of the required flow 
area and the dimensions of the connecting upstream component. 


Care should be taken to ascertain whether these assumptions 
apply for specific engine configurations. For example, a thin- 
walled cylinder subjected to an external collapsing pressure will 
fail at a much lower pressure than it would if it were subjected to 
an internal bursting pressure, as assumed in the duct weight cal- 
culation. If both ID and OD walls of the duct are exposed to 
ambient pressure, the ID wall should be sized to avoid collapse, 
such as determined experimentally by Stewart (ref. 2) for lap- 
welded steel tubes: 
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max 



1 t 

[i - > 

Jl-1600(§) _ 


or 


expressed in terras of rainimura wall thicknesaj 


( 30 ) 


^min 


D 

40 



1 - '^max 


689.5 


(31) 


WATE-S does not perform the above calculation to determine whether 
collapsing pressure sizes the ID wall. 


Crossover ducts are requited to join two centrifugal compres- 
sors, and are calculated in a different manner. The crossover 
duct diameter is assumed to be 1.6 times the exit tip radius of the 
upstream centrifugal compressor. The duct length is assumed to be 
20 percent of the diameter and the weight is estimated from the 
following expression: 


W 


crossover duct 


2.283 


(32) 


where R. is the exit tip diameter of the upstream centrifugal com' 
pressor. 


3.1,7 Burners 


The methods used for burners are based on a calculated volume 
of materials, similar to the previously described methods for 
ducts. However, wall liners, fuel manifolds, and fuel nozzle 
weights are also calculated as shown in Figure 9. The method 
described is used for primary burners, both axial and reverse- 
flow, as well as duct heaters and afterburners. 


The required Inputs include; 


Burner residence time 
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o Through-flow velocity 


o Burner mean radius for axial primary burners and duct 
heaters, and inner radius for reverse-flow primary bur- 
ners. 



Figure 9. Burner Schematic* 


The differences in burner types and configuration are 
reflected in the input values. When a primary burner is speci- 
fied, a frame weight may be added. Primary burners and duct 
heaters require a radius input, while the afterburner is assumed 
to have no inner wall. 


Burner flow area is determined from the input through flow 
velocity, the mean radius, and the inlet-corrected airflow from 
the engine cycle data. Burner length is found to give the speci- 
fied residence time based on the input velocity and entry condi- 
tions. Flow area is used to obtain the inner and outer dimensions 
of the burner (Rj^ and R^) from the specified mean radius. Inner 
and outer case thicknesses are determined from Equation (28). The 
assumed material is steel with an allowable stress of 48.3 KN/cm^* The 
we'ght of the inner and outer case is determined from Equa- 
tion (29) using the burner length L. 
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Liner weight la determined In a aimllar manner, assuming 
0,140"Cm thick steel walls, located at 20 percent of passage 
height from the Inner and outer case. The burner dome, fuel mani- 
fold, fuel nozzles, and other components are estimated by the fol- 
lowing equation developed by Boeing. 


"done - 

Total burner weight is the sum of the inner and outer cases and 
liners, burner dome and fuel nozzle system, and frame when speci- 
fied. 


3.1.8 Shafts 


A shaft is assumed to be the power connection component 
between compressors and turbines, or propellers/propfans and tur- 
bines. Figure 10 illustrates the shaft geometry and nomenclature. 


c 


Li 

I*-”'- L2 

^C ^ 

a 


n 



1 — 




ENLARGED SECTION OF CONCENTRIC SHAFT 
Figure 10. Shaft Schematic. 
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Multiple concentric shafts can also be specified, and will be 
sized around the inner shaft with OtSSl-cm radial clearance. 
Dimensions of the inner shaft are determined to provide the neces- 
sary torque capability at the specified allowable stress. Total 
shaft power Is the summation of work for all turbines on the 
shaft. Torque 1s calculated byt 

10 ^ 

T« — .-JCPW (34) 

iO 

where u) is the shaft rotational speed. 


Shear stress due to the torque load is defined by (ref. 2) 


T 


16 T Dq 

7t(D^^ - D 
o 



( 35 ) 


or in terms of the input diameter ratio (Dj^/D^) 




A gimllar procedure la uaed for concentric ehafts. The outer 
ahaCt's Inner diameter la found by adding 0.762 cm to Oq, and 
Equation (37) la aolved by iteration to aatlafy the dealred allow- 
able atreaa. 

While it is assumed In the shaft weight eatlmate that torque 
determines the shaft dimensions, It should be recognised that 
other design considerations may dictate shaft dimensions. Shaft 
critical speeds or longitudinal stiffness may actually design the 
shaft, but this is a function of bearing arrangement, mount stiff- 
ness, location of and stiffness of rotating masses. The calcu- 
lated shaft weight should be considered to be an absolute minimum, 
and can possibly be much larger when these other criteria are con- 
sidered. The rotational inertia of the shaft is not calculated, 
since It is a negligible quantity compared to the compressors and 
turbines. 

3.1.9 Frames 

A structural frame is normally required to span the engine 
flow path from the outer engine case to the shaft, usually to sup- 
port a bearing (as shown in Figure 11 for several typical 
engines). Mechanical design of the frame would require a defini- 
tion of all loads imposed on the frame under normal operating con- 
ditions, transients, and adverse operating conditions, such as a 
hard landing. This level of detail Is normally not available at 
the preliminary design stage for which WATE-S has been developed. 


Boeing had found, however, that the frame weight correlates 
well with the total frame-projected area. This data Is shown in 
Figure 12 for five types of frames commonly used: single-bearing 
frames with and without power takeoff (PTO) , turbine exit, and 
intermediate. Frame weight Is determined from this data, based on 
the local diameter and the type of frame specified. 
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Figure 12. Frame Weight. 

3.1.10 Nozzles and Thrust Reversers 

Unlike the rotating components, the loads and load paths of 
nozzles (particularly variable area C-D nozzles) are not readily 
defined on a general basis. A selected type of nozzle could be 
subjected to a detailed weight-estimating procedure, however, the 
trade-offs of internal and external performance with nozzle length 
and diameter would also be necessary to optimize the design. This 
type of data is not likely to be available at the level of develop- 
ment for which WATE-S is intended. 

A procedure has been developed that shows proper trends for 
multiple-stream nozzles and for variable geometry and fixed- 
geometry nozzles. Nozzle length is specified and should be 
selected to be representative for the type of nozzle. An effec- 
tive surface area is calculated based on the diameter of the con- 
necting component and the specified length. Only circular, coni- 
cal nozzles are assumed, and coannular nozzles can be represented 
by specifying a circular nozzle for each flow path. Plug nozzles 
can be represented by specifying a larger effective length; e.g., 
from nozzle entry to end of plug. 
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Wall tiilckrtfiSB la aisQumed to be a Q.l626"om composite of 
titanium ^^el3•-cm) and atalnleaB ateel (0.0813-om) below 644"K 
(1160*K) and 0.0813»om Btalnloaa Bteel above 644*K (1160”R}. 
Variable no^tBlea are calculated In the Bame manner except that the 
effective wall thlckneae la 2.75 times that of the fixed noseele. 

A weight estimating method for thrust reversers previously 
developed for aircraft preliminary design studies by Boeing was 
based on the weight of 18 different thrust reversers that are In 
current use. It had been found in correlations of this data that 

reverser weight (W) is a function of corrected mass flow (~^“) and 
nozzle pressure ratio ^nd is dependent on whether the stream 
is hot (primary) or cold (fan). The following relationship Is 
included In WATE-S: 


Hot (primary) 
C 


«1 

•^2 


^4 

0.52631 

2.2222 

191.87 

5.0348 

1.0036 

0.23014 

-0.5054 

0.56091 


The WATE-S method will apply the cold stream equation to a fan 
stream whether or not It Is heated by a duct burner. The hot 
stream equation is used for turbine exit or mixed-flow exhaust 
streams. 

3.1.11 Mixers 


A mixer is a mechanical device placed at the confluence of 
two coannular streams to Increase the mixing boundary so that 
thermal mixing takes place In a minimum length. Figure 13 illus- 
trates a typical example of a m'xer. This type of mixer is known 
as a daisy mixer or forced mixer. 
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figure 13. Daisy or Forced Mixer 


The required inputs include: 

o Mixer specific length, Lgp 

o Number of lobes or passages, N 

o Engine core and bypass flow areas. All and AIO 

Engine core and bypass annular flow areas are taken from the 
engine cycle data, and the inlet radius (R^) of the upstream com- 
ponent is used as a starting point for locating and R^, as shown 
in Figure 14. R^ normally will be the hub radius of the final 

turbine stage. 


L 

S 

yf 4A/m- 
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is based on core flow area (All) and 

A e (AlO + AII)/2 where AIO and All are inputs 

L = L/yfih/n is an input, calculate L 

sp 

N = Is the number of lobes or passages input 

Wmix “ 1-25 N (R^-R^) L 

ttR. ^ 

”plug “ SttR^L +sin' 15* ^2 

and K 2 are products of material density and thickness: 

= (0.127 cm) (0.00775 kg/cm^) = 0.00098 kg/cm^ (steel) 

K 2 “ (0.0813 cm) (0.00775 kg/cra^) = 0.00063 kg/cm^ (steel) 

Figure 14. Mixer Schematic. 

The mixer is assumed to consist of two components, the mech- 
anical mixer section and the mixer core plug, both fabricated from 
stainless steel with a density of 0.00775 kg/cm^. The average flow 




area of the mixer le used with the Input apeclflo length to deter- 
mine the length of the meohanloal mixer eeotlon. The meohanloal 
mixer eeotlon la aaaumed to be oonatructed of 0.121-cm material 
and the weight estimated from the following expreaaloni 


w 


mix 


. 1.25 N (Rq - R|„) bp (0.127) 


(40) 


The mixer core plug U assumed to be oonstruoted of 0.0013-om 
material and has, a oyllndrloal section as well as a conical sec- 
tion. The mixer plug weight la estimated from the following 


express! ton; 


^plug 


|2 7rRj^L + 



Sin 15** 


,j(0.08l3) 


(41) 


The total mix'*, weight is then: 


W 


w, 


mix 


+ W 


plug 


(42) 


3,1.12 Annulus Inverti ng Valve (AIV). 


This device has been used In some variable-cycle engines to 
invert the annular position of two concentric flow paths. It 
accomplishes the flow Inversion within a constant diameter enve- 
lope, and with constant-area duct passages. Figure 1 s ows a 
typical example of an AIV. This AIV was designed to vary the 
bypass ratio in a JT8D engine. 


The AIV weight method assumes a constructon similar to that 

shown in Figure 15, except that instead of ^ellw 

U1 is assumed to be a titanium honeycomb at 5.37 kg/m bel 

644®K (1160“R) and steel honeycomb at 9.13 kg/m above 6 

(1160'’R). If desired, different materials can be specified. 
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An empirical relationship, similar to the mixer method, had 
been developed by Boeing for estimation of the AIV weights 




where is the hub radius of the upstream connecting component, 
and and Rq are found to satisfy the input Mach number with the 
inlet corrected airflow (see Figure 16). The number of passages 
(N) is an input, and material weight per unit area (W/A) is 
selected depending on the stream temperature. Length (L) of the 
AIV is calculated from the input specific length, Lgpt 

L = Lgp/4A/7 (44) 


^OUTER I— ^ 


^INNER 

R. I 



OUTER CYL 


INNER CYL 


Lsp - - INPUT A - 

FROM MACH INPUT 

^INNER ^UTER FOUND 

L IS DETERMINED 
WTIC-2 R, L WTSI 
WTOC-2 Rq L WTSO 

WTWAUl-IK, Rm + >^2 N(Ro-R|h ^ WTSW 

WTSI. WTSO AND WTSW ARE MATERIAL WEIGHT PER UNIT AREA 

K^- 3.927 K2-1.26 

Figure 16. Annulus-Inverting Valve Schematic. 





Specific length la preferred «a an Input because It la non- 
dlmenalonal, and It la a major variable that determlnea AIV pres- 
sure loss. A relatively good compromise between size and perform- 
ance la achieved when N"8 and Lgp*0.8 to 1.0» which results in a 
pressure loss between 2.5 and 1.5 percent. 

If the AIV is of the switching type, where one half Indexes 
In a rotational direction relative to the other half to change 
flow-path orientation, an actuator weight Is estimated at 10 per- 
cent of total AIV weight. Additional structure to support the 
rotating n^lf Is not Included and should be represented as an 
additional frame. 

3.1.13 Gearboxes and Transmissions 

A method of estimating the weight of various types of gear 
systems has been previously developed by Schmidt (ref. 3). This 
method has proved applicable for estimating the weight of turbofan 
engine gearboxes, such as the one in the Garrett TPE731. The 
equation used for total gearbox weight Is 

Gearbox Wt » 15.5 . 

„here maximum torque transmitted 

and GR = gear ratio (> 1.) 

The gearbox used In turboprop engines la weighed with a dif- 
ferent method. 

Gearbox Wt • 471.86 
3.1.14 Heat Exchangers 

Both rotary and fixed heat exchanger weights can be esti- 
mated. Methods previously developed produce adequate results for 
preliminary design purposes, see Figure 17. 
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ROTAKV HIAT CKCHANOER 


Figure 17. Heat Exchangers. 


For rotary heat exchangers, a ceramic core is assumed. 
Weights of thJs type of heat exchanger have been determined by the 
corning Glass Company (ref. 4) and are represented in Table 1 for 
various levels of effectiveness and pressure loss. This data is 
developed for a total corrected airflow of 90.7 kg/s. For other 
sizes, these weights are scaled directly with corrected flow. 


TABLE I. CERAMIC ROTARY REGENERATOR WEIGHT 


90.7 kg/s corrected weight flow 


BPR 

3 

7 


iP/P {%) 

5 10 

5 10 

hBdB 

Weight (kg) 

e a 80% 

85% 

90% 

305.7 245.8 

440.4 342.0 

735.7 538.9 

284.9 201.8 

414.6 290.3 

647.7 492.1 

272.2 191.9 

378.3 279.9 
595.6 450.9 




Pixer^-tube heat exchangers are estimated by a heat-'tranaf er 
analysis (ref. 5) where the required tube surface area is found to 
give the specified effectiveness. Plow area of the tubes is found 
from an input Mach number# number of tubes, and corrected flow. 
An estimate of the pressure drops Is made for the input number of 
tubes based on input relative roughness values for each flow path. 
The flow Reynolds number is calculated and the Moody diagram is 
used to estimate the friction factor, and thus the pressure drop. 
These calculated AP/P's can be used to ascertain whether the num- 
ber of tubes selected is reasonable. Wall thickness of the tubes 
is determined by Equation (28) to satisfy an assumed allowable 
stress of 34.5 KN/cm^ and a density of 0.00465 Kg/cm^ below 644®K 
(1160®R). A stress of 48.3 KN/cm^ and a density of 0.00792 Kg/cra 
is assumed above 644®K (1160®R). Minimum wall thickness of 

0.0254-cm is used for the tubes. The length of the tubes is deter- 
mined to satisfy the surface area requirements. Fixed-tube heat 
exchanger tube weight , is then found by 

“tubes = 

where and R^^ are the tube radii and L is the total length 
required. A wrap-up factor of 1.85 is used to account for the 
weight of the casings, mounting hardware, manifolds, and other 
equipment that may be necessary. 

^HTEX “ * '^tubes 

3.1.15 Prop ellers and Propfans 


The following expressions for estimating the weight and diam- 
eter of advanced technology propellers and propfans were derived 
from methods developed by Hamilton Standard (ref. 6). It is noted 
that these expressions are based upon a conventional blade of the 
stated materials with a generally rectangular platform shape, and 
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a atandfi.:i Uemgaa shank. The blade weight relatione are for com- 
mon propeller materials? composite and aluminum. If other 
materials are desired, a density ratio can be used to estimate the 
weight. However, allowable working stresses and blade stiffness 
may vary and considerable judgement is required to account for 
these effects. 

Thp expression used in WATE-S for double acting, advanced 
technology propellers is 



/ n / 

M /a 

F / Ut ) 

Prop Wt = K 

^17048j ( 

V 

00 J \319.2/ 


,. 0.5 / PW \°*^^1 

(Mn +1) I “5 I 1 

\° / 


where K =42.6 for composite materials 
57.6 for aluminum 
D = diameter (m) 

N * number of blades 
AF » blade activity factor 
PW = maximum shaft power (KWatt) 

Mn » design cruise Mach number 

= 100-percent design tip speed (m/s) 


The expression for propfans is 



The input quantities are tip speed, PW/D ‘ or D, the number of 
blades, activity factor, and design cruise Mach number. The cal- 
culation procedure determines the diameter if not input, then the 
propeller RPM is calculated from the tip speed, 
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RPM 


60 
7T D 


(51) 


and finally the propeller or propfan weight Is calculated from the 
above expressions. 

3.1.16 AcceaBories 


Accessory weight data was collected from twelve general avi- 
ation aircraft manufacturers and was supplemf 'ted with Garrett 
data for several engines. The data Is presented In Table IX and 
includes 22 various turbofan, turboprop, and turboahaCt engines 
ranging In sea level static thrusts from 2.67 KN to 22.46 KN and 
maximum shaft horsepowers from 0.373 mWatt to 1.163 mWatt. 

A tremendous amount of scatter exists in the data, and hence, 
separate correlations for each accessory were not possible. Bat- 
tery weights varied and were generally the heaviest accessory. 
For these reasons, the starter/generator , hydraulic pump, oil 
pump, and fuel pump were lumped together and a single correlation 
derived for these accessories (note that the battery Is not 
included) , 

ACCS Wt = 18.2 + 0.02 (bare engine weight) (52) 

If this calculation results in an accessory weight fraction less 
than 10 percent of the bare engine weight, then 10 percent of the 
bare engine weight is used for the accessories. The bare engine 
weight where this break point occurs Is 198.4 kg. 

3 . 2 Other Program Functions and Capabilities 

3.2.1 Flight Envelope Maximization 

In the normal use of the WATE-5? program, a flight envelope of 
engine cycle data will be generated. Since the weight of each 
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•Op* battery for A/C 

••Double-actiPQ/Ei-'^l^-acting prop. 







































































aomponent la affected by its maximum wock, flow, temperature, and 
apaedj these maximum values are stored for use in the dimension 
and weight calculations. The flight condition is given in the 
output data where the maximum condition occurs for each component. 

3.2.2 Design Limits 

As an aid to assist the user In achieving a reasonable engine 
design, the output will provide a warning and suggested corrective 
action to bring the engine design within reasonable limits. These 
limits can be specified, or default values will be used if not 
specified. Table III shows a list of warnings and corrective 
actions , 

3.2.3 Automatic Airflow Scaling 

The WATERS program will automatically scale the engine ±20 per- 
cent of the size that is defined by the thermodynamic input. 
Up to six selected scale factors can also be specified. A scaling 
exponent (e) is calculated for each scaled engine based on the 
following expression: 



The scaling exponent (e ) for each engine size is provided in the 
output data. 

3.2.4 Engine Center of Gravity 

The center of gravity of each component, except the rotating 
components, is assumed to be the midpoint of its length. The 
moment relative to the front flange of the engine is determined as 
a function of the position of each component. 
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Warning Meaaago 

Action Recommended 

Default Value 
Tested Against 

Blade Centrifugal 
Stress Exceeded 

Reduce shaft speed 
(by rpm scalar) or 
increase exit Mach 
number 

fr m 34.5 KN/cm? psi HPT 
(f • 41.4 KN/cm, pal LPT 
<r m 55,2 KN/om pal 
fans and axial 
compressors 

h/t TOO Large 

Reduce hub/tip 
ratio input 

h/t » 0.93 HPC exit 

h/t Too Small 

Increase hub/tip 
ratio input 

h/t - 0.32 fan, LPC, HPC 
inlet 



h/t - 0.50 HPT, LPT exit 

— 
Turbine Work 
Too High 

Add turbine stages 
decrease input 

. ^H » 139.6 KJ/kg 
e HPT , LPT 

Stage PR Too 
High 

Reduce stage pres- 
sure ratio input 

PR ■ 1.8 fan 
PR « 1.65 HPC, LPC 

Flow Velocity 
Too High 

Decrease stage 
inlet Mach number 
input 

Mn = 0.60 all 

components 

Blade Size Too 
Small 

Change overall 
pressure ratio or 
reduce h/t input 

h„ ■ 1.016 cm all axial 
“ compressors 

Compressor Work 
Coefficient Too 
High 

Add compressor 
stages or increase 
hub/tip ratio 

= 0.9 all axial 
compressors 





For rotating components, the weight of each stage Is ass’mo^ 
to act at mid'-length of the stage t Moments are summed about the 
front flange for each stage. The center of gravity of each rotat" 
Ing component and the total engine 1s calculated. CG locations 
are shown In the output data for the total engine and each com'' 
ponent . 

3 . 3 Program Validatio n 

A verification of the accuracy of the WATE-S method can only 
be done by applying It to various types of engines and comparing 
the results with the actual measured engine weight and dimensions, 
or with those estimated by the manufacturer for proposed engines. 
Since the manufacturer's estimate of proposed engines also 
includes some error, the real deviation or error of the WATB-S 
method can only be found by comparing engines that have been built 
in production quantities. 

In order to demonstrate the accuracy of WATE-S, the NASA Pro- 
gram Manager and Garrett jointly selected four Garrett propulsion 
engines for comparison. These included both production and pro- 
posed engines. Results of this comparison are presented in 
Figure 18. As can be seen, the predicted engine weight and over- 
all dimensions of the selected engines are within the ±10 percent 
accuracy goal. 




4.0 USKKS MANUAL 


This section contains a description of the WATE-S input and 
output data, values of typical inputs, sample cases, and program 
structure definition. WATE-S is designed to function with a com- 
ponent type thermodynamic engine cycle analysis program. The 
thermodynamic engine design point may be used to generate the 
WATE-S inputs, or additional off-design points can be used and the 
maximum conditions of component work, flow, temperature, and pres- 
sure are used to size the turbomachinery components. In order to 
achieve the most accurate engine weight estimate, the off-design 
cases should include the maximum performance levels required of 
each component. 

WATE-S will also accept input weight scalars for each com- 
ponent so that components may be selectively eliminated, 
increased, or dectca&ed in v;clght to determine sensitivities, etc. 
Input and output units may be either English or SI. 

4.1 Input Description and Formats 

The inputs of WATE-S ate arranged in three NAMELIST groups: 
$C for the configuration data, $T for the thermodynamic data, and 
$W for the aeromechanical design data. The configuration data is 
input only once per WATE-S execution, whereas the thermodynamic 
and aeromechanical design inputs are grouped in pairs. Table IV 
provides a summary of the input variables. Figure 19 illustrates 
a typical card stacking arrangement to execute WATE-S, The 
following sections discuss each of these NAMELIST groups in detail 
and provide typical and default values where appropriate. 
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4 . X . 1 Configuration Inputa 


The engine configuration input is contained in the $C name" 
list. It is preceded by a title card describing the engine con- 
figuration or any other description the user might desire. The 
configuration inputs fully describe the component connectivity and 
define the flow stations throughout the engine. The following 
variables were originally supplied by NNEP (ref. 7) in WATE-2, but 
in WATE-S are required input. 


"NCOMP" is an integer variable defining the total number of 
components used in the WATE-S engine simulation. 

"NOSTAT" is an integer variable defining the total number of 
flow stations used for the thermodynamic simula- 
tion of the engine. 

"JTYPE" is a one-dimensional integer array containing the 
type of component. 


1 

Inlet 

•INLT' 

2 

Duct 

•DUCT' 

3 

Water Injection 

'WINJ' 

4 

Compressor 

•COMP' 

5 

Turbine 

•TURB' 

6 

Heat Exchanger 

'HTEX' 

7 

Splitter 

•SPLT' 

8 

Mixer 

'MIXR' 

9 

Nozzle 

'NOZZ' 

10 

Propeller 

'PROP' 

11 

Shaft 

•SHPT' 


(this component known as 'LOAD 
in NNEP) 
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"JCONF" la a two-dlroenalonal integer array defining the 
component connectivity ualng the flow atatlon 
number a. It la of the form JCONF(M,N) where M la 
the component number ^ and N la the variable number 
aa deacrlbed below. 


Location Description 

1 Is the primary upstream airflow station number 

for flow components) or the first component 
hooked onto a shaft, or the shaft hooked onto 
a propeller. 


2 is the secondary upstream station number, or 
the second component hooked onto a shaft. 

3 is the primary downstream station number, or 
the third component hooked onto a shaft. 

4 is the secondary downstream station number, or 
the fourth component hooked onto a shaft. 


A default engine configuration has been included in WATE-S. 
Figure 20 illustrates this engine; a two-spool turbofan with axial 
single-stage fan, axi-centrif ugal HP compressor (three-stage 
axial), axial single-stage HP turbine, and axial two-stage LPT. 
The default $C namelist data corresponding to this engine config- 
uration are shown in Table V. 




1 


SB 



i 



FAM 


OKtGINAL PAGE IS 
OF POOR QUALITY 
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Figure 20. Default Engine Configuration 


























The thermodynamic input is contained in the $T namelist. It 
Is preceded by a title card describing the engine being weighed or 
any other description the user might desire. The thermodynamic 
input fully describe the properties for each flow station identi- 
fied in the configuration input. These properties include mass 
flow rate, total pressure, total temperature, and fuel/air ratio. 
These Inputs are used as component upstream flow conditions. 


is a one-dimensional array containing the physical 
mass flow rate at each station (Ibm/s) . 


"TOPRES" is a one-dimensional array containing the total 
pressure at each station (psi) • 

"TOTEMP" is a one-dimensional array containing the total 
temperature at each station (*R) . 

“PAR" is a one -dimensional array containing the fuel/air 

ratio at each station. 


in addition to the above, three additional arrays of data are 
required. These arrays were originally filled by the NNEP code in 
WATE-2, however, since WATE-S is not configured with an engine 
cycle code, separate input variables were needed. 

"PERPP" is a one-dimensional array containing twelve 
values which describe general performance param- 
eters for the engine. 
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6 not thruat/lnlot total flow 

7 total inlet drag 

8 total brake abaft power (bp) 

9 inatalled tbruat 

10 installed TSFC 

11 spillage and lip drag 

12 boattail drag 

Tbia data is used for information purposes only, and bence, is 
considered optional in WATE-S. 

"DATINP" is a two-dimensional array containing several mis- 
cellaneous component properties. It is of the form 
DATINP (N,M) where M is the component number, and N 
is the variable number as described in Table VI. 
The only required inputs in this array are the 
shaft gear ratios. The remainder of the inputs are 
optional. The bleed fractions will be calculated 
from the "WTP" array input and the Mach number and 
altitude are for information purposes only. 

"DATOUT" is a two-dimensional array containing several 
other miscellaneous component properties. It is of 
the form DATOUT (N,H) where M is the component 
number, and N is the variable number as described 
in Table VII. The required Inputs in this array 
Include the turbine rpm, if the turbine is not con- 
nected to a compressor (free-turblne turboprop and 
turboshaft engines) , splitter bypass ratio, and 
mixer primary and secondary flow areas. The 
remainder of the array is calculated from the 
"WTP”, "TOPRES", TOTEMP", and "FAR" arrays, if not 
input . 
















































4,1,3 Aegomechanloal Input 


The aegomechanloal Input Is contained In the $W namelist. 
These variables control the aerodynamic and mechanical suing oC 
the components. The following sections describe the aero- 
mechanical input and give the default values where applicable. 

4. 1.3,1 Miscellaneous Variables 

"ACCARM" is a variable that contains the value of the cen- 
troid distance for the accessories component in the CG calcula- 
tions. If no value is input r accessories are not included In 
center-of-gravity calculations. 

"DISKWI" is a variable that is used as an indicator for the 
axial compressor disk weight method. (Default ^ i.) 

0 - Do disk weight calculations using the Boeing WATE-2 

method. 

1 - Do disk weight calculations using the Garrett 

WATE-S method. 

"ACCS" is a one-dimensional array containing three values 
which control the accessory weight calculations. (Default >40., 


0.02, 0.10) 



ACCS(l) 

a 

Fixed accessory weight (Ibm) 

ACCS (2) 

a 

Accessory weight fraction of bare engine 
weight 

ACCS (3) 

a 

Minimum total accessory weight fraction 
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"ISCAEtE" ia a one-dimenolonal intogor array containing thraa 
vaXuoe which control tha engine acallng Xogio ot the program. 
(Default " 1, 3, X) 


ISCALE(X) 


Output indicator for scaled engines. 
Unsealed engine output set by "lOUTCD". 


0 - Short form engine weight, length, and maximum 

radius. 

1 - Long form component weights and dimensions. 

2 - Debug option. 


ISCALE(2) Number of scaling points. 

ISCALE(3) Not used. 


"SCALE" is a one-dimenaionaX array containing six values. 
The values correspond to the engine scaling factors desired. The 
number of "SCALE" values must be equal to the value of ISCALE(2). 
First value must always equal X. (Default » 1, 0.8, 1.2, 3*0.) 

"IWT" is a integer variable to control the execution of the 
component weight calculations. (Default ■ 1) 

0 - Do not do weight calculations. 

1 - Do weight calculations without the thermodynamic 

parameter maximization feature of the WATE-S cod^. 
This is used for design point paramteric weight 
analysis . 

2 “ Do weight calculations using maximum thermodynamic 

parameter maximization feature. 
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3 - DO weight oaloulfltion with airflow acallng without 

the thermodynamic parameter maximlssation feature, 

4 - Do weight calculation with airflow acallng using 

the thermodynamic parameter maxlraization feature, 

"IPLT" is a logical variable to control the gas path printer 
plot option of WATE-S. (Default - T) 

T - Gas path layout 

p - No gas path layout 

”ISir* is a logical variable to specify the units of the $W 
namelist variables. (Default “ P) 

T - SI units input 

P - English units input 

"ISIO" is a logical variable to specify the units of the 
output of the weight calculations. (Default » F) 

T - SI units output 

F - English units 

"lOUTCD” is an integer variable to control the amount of 
output generated by the component and engine weight calculations. 
(Default • 2) 

0 - Short form engine weight, length, and maximum 

radius. ' 

1 ~ Long fori., component weights and dimensions. 


2 


Debug option. 



4 • 3 • 2 Length Contributing Veofeor 

"ILRNG" is a one-dlmonslonal Integer array specifying only 
those components which contribute to the total additive engine 
length. The component numbers are specified in the order that the 
components would add In length to achieve the total engine length. 
This must start with the inlet and end with the furthest down- 
stream nozzle. The default "ILENG” values for the default engine 
configuration illustrated in Figure 20 aret 

ILENG = 1, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 

The default "iLENG" input does not include the bypass flow compo- 
nents {COMP 2, DUCT 3, NOZ 4), reverse-flow burner (DUCT 10), or 

shafts (SHFT 16, SHPT 17) since these components do not contribute 
to the total engine length. 

4. 1.3. 3 Mechanical Design Indlcahrtrg 


"IWMEC is a two-dimensional integer array containing the 
mechanical design indicators, it is of the form IWMEC(N,M), where 
M is the component number, and N is the variable number as 
described in the following sections for each component type. The 
mechanical design indicators must be specified for each component. 
Table VIII illustrates the default "IWMEC" values, which are 
representative of the default engine (see Figure 20) . 








































Location 

1 


2 


3 


I CompreasotB 


Peace Iptlon 


Type of corapEeaaor being weighed 


(1) 

•FAN' 

- Typical fan 




(2) 

•FO' 

- Outer portion 
ter fan 

of 

nonrotating split- 

(3) 

•FI' 

- Inner portion 
ter fan 

of 

nonrotating split- 

(4) 

•RSFO' 

- Outer 
fan 

portion 

of 

rotating 

splitter 

(5) 

•RSFI' 

- Inner 
fan 

portion 

of 

rotating 

splitter 

(6) 

•lpC 

- Low-pressure compressor 


(7) 

•HPC 

- High- 

pressure 

compressor 



Thia indicates if the compressor has stators or if 
the compressor is a centrifugal compressor. 

0 - Stator weight is not calculated 

1 - Stator weight is calculated 

2 - Centrifugal compressor 

This is the indicator for 'front' frames in com- 
pressors. This input may be: 


0 - No frame 

1 - Single bearing frame without Power 

Takeoff (PTO) 

2 - Single bearing frame with PTO 

4 - Two bearing frame, such as the frame 

in front of the HPC in the JT8D tir 
JT9D, which extends outward to the 
fan outer case and holds two bearings 
with PTO 



Deacr Iptlo n 

This is the Indicator for the 'rear' frame In a 
eompreaaor 

0 “No frame 

1 “ Single bearing frame without Power 
Takeoff (PTO) 

2 “ Single bearing frame with PTO 

4 “ Two bearing frame, such as the frame 

in front of the HPC in the JT8D or 
JT9D, which extends outward to the 
fan outer case and holds two bearings 
with PTO 

5 This is the component number connecting to this 
component for split flow compressors only. If 
this is the Pan Outer, the Pan Inner must be 
specified. If this is the Rotating Splitter 
Outer, the Inner Splitter must be specified, and 
vice versa. If {-) , a dummy compressor is 
assumed. 

6 Gearbox indicator. Input is component number of shaft, 
0 for no gearbox. 

7 Number of stages; if 0, will calculate number of stages 
assuming equal work based on the maximum first-stage 
pressure ratio input (DESVAL(2,M)]. 

4. 1.3. 3. 2 Turbines 


Location 

4 


Location 

1 

( 8 ) 

(9) 

2 


Description 

This is the type of turbine 

'HPT' - High-pressure turbine 

'LPT' - Low-pressure turbine 

Indicator for turbine exit frame 
0 - No frame 

3 ~ Single bearing turbine exit frame 

for large engines 

S '• Small turbine exit frame 
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LQoatton 

3 


4 


5 


6 


7 


Deacrlptlon 

Compresaoi: number from which the rpm is deter- 
mined; i£ 0, DAT0UT(2,M) is used for turbine rpm 
(axial turbines only) • 

Component number from which the mean radius limit 
for the turbine is determined. If the component 
number is positive, the outlet outer dimension is 
used. If negative, the inlet outer dimension is 
used. If 0, it will use the outlet of the feeding 
component. 

Number of stages; if 0, will calculate number of 
stages assuming equal work based on turbine mean 
radius limit and mean work coefficient. 

Indicator for axial or radial turbine 

0 - Axial turbine 

2 ~ Radial turbine 


Not used. 


4. 1.3. 3. 3 Burners 


Location 

1 


2 


Description 



This is 

the type of burner being 

weighed 

(10) 

'PBUR' 

- Primary burner 


(11) 

'DBUR' 

- Duct burner (a mean 
fled) 

radius is speci- 

(12) 

'AUG' 

- Augmentor (no inner 

wall) 


This is 
only for 
bearing . 

the Indicator for frame 
primary burners. This 

weight, normally 
frame includes a 


0 

- No frame 



1 

- Frame 



3-7 


Not used. 


4. 1.3, 3. 4 PuctB 
Location 

1 

2 


PeBcriptlon 

(20) 'DUCT' 

Indicator as to type of duct 


1 

2 

3 


- Dummy - i.e.i no weight or length 

- Length input 

- Length derived as in a duct connect- 
ing a splitter and a mixer 

- Crossover duct for centrifugal com- 
pressors 


3-7 

4. 1.3. 3. 5 Shafts 
Location 


Not used. 


Description 


1 (13) 

2 

3-7 

4. 1.3. 3. 6 Mixers 
Location 
1 


•SHAF' 

Shaft number from inner to outer. 
Not used. 


Description 


(14) 

(15) 


Type of mixer 
•MIX' 


•PMIX 


- The coannular emergence of two 
streams without mechanical mixer 

- Forced mixer^ mechanicalf i.e.» Daisy 
lobed mixer 


3-7 


Indicator for primary input node 

0 Engine core is inner (primary) 

1 Engine core is outer (secondary) 

Not used. 
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4.1.3.3.7 Nozzles 
Location 


Pescr lotion 


1 

2 


5-7 


(16) 'NOZ' 

Nozzle type 

1 - Convergent 

2 - C-D variable area 

Component number from which the nozzle inlet dlam 
’iVL “"f K this dlliete? i,‘ 

m}nL a. ^ ^ P^^®vious componont deter- 

mines the diameter, this location may be zero. 

Thrust reverser type 

0 - None 

1 - Pan 

2 - Primary 

Not used. 


4. 1.3. 3.8 Heat Exchangers 
Location 


Description 


1 

2 


(19) 


•HTEX' 

Heat exchanger type 

1 - Fixed tube 

2 - Rotary 
Plow Direction 


- Parallel flow 

- Counter flow 


4-7 


Not used. 


4. 1.3. 3. 9 

Solitters 



Location 



Description 


1 

(17) 

' SPLT ' 



2 


0 

- Inner stream is primary 




1 

- Inner stream is not primary 


3-7 


Not us3d. 



4.1.3.3.10 

Annulus Inver tine Valve 

1 

1 

Location 



Description 


1 

(18) 


'VALV 


2 


Location 

of Valve 




1 

- Inner 




2 

- Outer 


3 


Component 

number of opposite duct 


4 


0 if fi::ed, 1 if movable 


5-7 


Not used. 



4.1.3.3.11 

Propellers and PronfanR 


Location 



Description 


1 

(21) 

•PROP’ 



2 


Propeller 

type 




0 

- Hamilton Standard double 
advanced technology aluminum 
ler . 

acting, 

propel- j 



1 

- Hamilton Standard double acting, 
advanced technology composite propel- 
ler. 



2 

- Hamilton Standard propfan. 

] 

3-7 


Not used. 


1 

1 
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Aeromechanlgal Design VarlahiAH 


"DESVAL" ia a two-dimensional array containing the aero- 
dynamic and mechanical design inputs. It is of the form 
DESVAL(N,M), where M is the component number, and N Is the vari- 
able number as described in the following sections for each com- 
ponent type. The default values are stored in a separate array 
(DEPAUL) which cannot be altered by program WATE-S input, a sum- 
mary of these arrays is contained in Table IX. The default values 
of the "DEPAUL" array are illustrated in Table x. The "dbpaul" 
array values are used when the "DESVAL" input is not specified for 
a given component. Table XI contains typical ranges of the 
"DESVAL" input variables for each component type. 



















































































TABLE X. DEFAUL ARRAY VALUES 
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4. 1.3, 4.1 CompgeBiBoc 


Arr ly 


Location 


Description 



1 


Compressor Inlet Mach number. 



2 


Maximum first-stage pressure ratio* 
If 1WMBC(7,M) equals 0. 

required 

3 


Compressor Inlet hub/tlp ratio. 



4 

Axial ; 

Blade tip solidity* ratio of blade tip chord and 
blade spacing. 


Centri- 

fugal: 

Number of blades. if {-) * splitter 
calculated. 

weight 

is 

5* 


Blade aspect ratio of first stage. 



6* 


Blade aspect ratio of last stage. 



7 


Compressor exit Mach number. 



8 


Maximum compressor inlet temperature 
design point temperature is to be 
material selection. *R (**K) 

. ZERO 
used 

if 

tor 

9 


Maximum compressor outlet temperature. ZERO 
if desired point temperature is to be used for 
material selection. ®R (®K) 

10 


Maximum speed ratio - . 



11* 


Blade material density. . ZERO i£ WATE-S is 
select material. Ibrn/in"^ (Kg/cm'*) 

to 

12 

Axial: 

Compressor design type 





1. Constant hub radius design. 





2. Constant mean radius design. 





3. Constant tip radius design. 




Centri- 
fugal : 





♦Not used for centrifugal compressors. 
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LooatlQn 

13 

14* 

15 

16* 

17* 


Deaorlptton 

RPM acaler, normal input la 1. - use to match 
known rpm or modify tip apeed correlatlona. 

Temperature at which a change of material la 
required. ®R (®K) 

Comproeeor weight scalar r Input ZERO if no 
scaling Is desired. 

Blade taper ratio. 

Blade volume factor. 


"^^ot used for centrifugal compressors. 
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4 . 1 « 3 . 4 . 2 Turblneo 
l^oc at ion 


Dr acription 

1 Turbine inlet Mach number. 

2 Axial: First-stage mean work coefficient - X ■ 

2 

U ^ 
m 

Radial: Radial turbine work coefficient - 


3. 


4* 

5* 

6 

7 


8 * 


9 

10 * 


Axial: Blade tip solidity, blade tip chord/blade 

spacing. 

Radial: Number of blades. if (-), splitter weight is 

calculated. 

Blade aspect ratio of first stage. 

Blade aspect ratio of last stage. 

Turbine exit Mach number. 

Disk reference ftress - 0.2-percent yield, 
lbf/in2 (Newton/ cm2) . 

Turbine design type 

1. Constant tip radius design. 

2. Constant mean radius design. 

3. Constant hub radius design. 

Maximum speed ratio - ^’*Mj^ax'^^^^design * 

Turbine control radius, inch (cm) - blank if 
transferred from a component. This overrides 
IWMEC input if nonzero. 


*Not used for radial turbines. 
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Location 

Description 

11 

Density of material in turbine blades, Ibm/ln^ 
(Kg/cmJ) , 

12 ^ 

Blade volume factor. 

13-14 

Not used. 

15 

Turbine weight scalar. Input ZERO, if no scal- 
ing is desired. 

16* 

Turbine blade taper ratio. 

17* 

Stator blade volume factor. 


*Not used for radlnl turbines. 


4. 1.3.4. 3 Burners 


Location 

1 

2 

3 


4 


5 


Description 

Burner through-flow velocity, ft/s (m/s) . 

Burner airflow residency time, s. 

Burner mean diameter, inch (cm). If zero, diameter 
is calculated to match a connecting component (see 
below) . 

Component number for calculating burner reference 
diameter, if other than upstream component. Not used 
v’hen diameter is specified. If greater than zero, 
connecting component mean diameter is used for burner 
mean diameter. If less than zero, connecting com- 
ponent outer diameter is used for burner inner diam- 
eter (reverse-flow burners) . 

Number of cans for can burners. 


6-14 Not used. 

15 Burner weight scalar. Enter ZERO, if no scaling ip 

desired. 


16-17 Not used. 
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4. 1.3. 4. 4 

L ocatio n 

1 

2 

3 

4 


Quota 


Deacri ption 

Duct inlet Mach number. 


Length to height ratio of 
IWMEC(2,M) eguala 2. 


duct, 


required if 


Duct mean diameter, inch (cm), if o, duct diameter 
Is cslculstsd to mstch d connectlnQ componsntt 

Component number for calculating mean duct diameter 
Not used when mean duct diameter is specified! 
Enter -l, if upstream connecting component is to be 

USSu « 


5-14 Not used. 

Weighu scalar. ZERO, if no scaling is desired. 
16-17 Not used. 


4.1. 3.4.5 Shafts 


Location 

1 

2 

3 

4-14 

15 

16-17 


Description 

Shaft allowable stress, Ibf/in^ (Newton/cm^) . 

Shaft material density, Ibm/in^ (Kg/cm3) . 

Diameter ratio of shaft, D, /D ^ . not used for 

outer shafts. inner^ ‘"outer ' 

Not used. 

Shaft weight scalar. ZERO if no scaling desired. 

Not used* 
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4.1*3«4t6 

Location 

1 


2 

3-14 

15 

16-17 


Mixers 


Description 


Effective ^ngth-to-dlemeter tetlo 
mlxet, LZ/vS, where L Is the mixer length Inlet to 
exit, A is the total flow area. Enter 0, if not a 
mechanical (forced) mixer. 


Number of passages (or lobes) in mixer. 

Not used. 

Weight scalar. Enter ZERO, if no scaling desired. 
Not used. 


4. 1.3. 4. 7 

Location 

1 

2 

3-14 

15 

16-17 


Nozzles 

Description 

Length to diameter ratio of nozzle. 

Bypass ratio for mixed flow nozzle for thrust 
reverser weight. 

Not used. 

Weight scalar. Zero, If no scaling desired. 

Not used. 


4. 1.3. 4. 8 

Location 

1 

2 

3 

4 

5 

6 

7-14 

15 

16-17 


Heat Exchangers 

Description 

Number of tubes if "Fixed" type. 

Mach number in primary stream (cold side). 
Mach number in secondary stream (hot side). 
Engine bypass ratio if "Rotary" type. 

E/D for primary stream if "Fixed" type. 

E/D for secondary stream if "Fixed" type. 
Not used. 

Weight scalar. 

Not used. 


Zero if no scaling desired. 



4 , 1 , 3 . 4 . 9 Split ters 


Location 


3-14 


D esc r iptl on 


Only Input 1£ £lrst calculated component In £low 
path. Inlet Mach number. 

Inlet hub/tip radius ratio. 

Not used. 


16-17 


Weight scalar. ZERO, i£ no scaling desired. 
Not used. 


4.1.3.4.10 Annulus Inverting Valve 


Location 


Description 


9-14 


16-17 


Specific length - ratio of length to effective diam- 
eter of the AIV, lV4A/7t 

Number of passages. 

Mach number of inner passage. 

Mach number of outer passage. 

Hub radius in Inches (cm) ; or component number from 
which hub radius is taken; or blank, if feeding com- 
ponent determines the hub radius. 

Inner cylinder weight - Ib/ft^, (Kg/m2) , 

Outer cylinder weight - Ib/ft , (Kg/m2) . 

Wall weight - Ib/ft^, (Kg/m2) . 

Not used. 

Weight scalar. ZERO, if no scaling desired. 

Not used. 




4. 1.3. 4. U 


Location 

1 

2 

3 

4 

5 


Propel lera and PcopCans 

Descript ion 

Design cruise Mach number of aircraft. 

Number of blades. 

Activity factor. 

Tip speed - ft/a (m/s) . 

Ratio of propeller shaft power to diameter squared 
hp/ft2, (KW/m^). 


6-14 Not used. 

15 Weight scalar. Zero, if no scaling desired. 

16-17 Not used. 


4 . 1 . 3 , 5 Design Limits 

"DESLIM" is a one-dimensional array containing the design 
limits for WATE-S. Table XII describes the design limits and the 
default values currently in the code. These can be changed as 
desired through the namelist $w input. If these limits are 
exceeded, the component weight and dimension calculations continue 
and a warning message is printed out. 

4.2 Pro gram Output 

The output of WATE-S is contained in three sections, coincid- 
ing with the input: configuration, thermodynamic, and aeromechani- 

cal. The configuration output contains the $C namelist Inputs and 
a simple configuration layout. The thermodynamic output contains 
the $T namelist inputs. 

The aoromechanlcal output may be selected in either English 
or SI units. The units used are presented in Table XIII. The type 
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TABLE XIII. OUTPUT UNITS . 


VARIABLE 

SI UNITS 

ENGLISH UNITS 

Velocity 

m/a 

ft/s 

Temperature 

®K 

®R 

Pressure 

N/m^ 

Ibf/ft^ 

Area 

m 

ft^ 

Stress 

N/cm^ 

Ibf/in^ 

Density 

Kg/cra^ 

Ibm/in^ 

Weight 

Kg 

Ibm 

Length 

cm 

in 

Enthalpy 

Joule/Kg 

Btu/lbm 

Power 

KWatt 

hp 

Weight flow 

Kg/s 

Ibm/s 

Weight flow/unit area 

Kg/m^s 

Ibm/ft'^s 

Radius 

cm 

in 



of miita In use are noted in the units section of the aero- 
mechanical output. The level of output is controlled by the $W 
namelist integer variable "lOUTCD" and has three output options. 
Examples of these output formats are shown in Section 4.3, Sample 
Cases for two different engines. 

A flow path printer plot is also available when $w namelist 
logical variable "IPLT" is true. This plot is automatically 
scaled to fit on one page of output. Examples of this option are 
also contained In Section 4.3, Sample Cases. 

4.3 Sample Cases 

Two sample cases are included herein to illustrate the 
execution and output of WATE-S. The first is a single-spool 
turboprop and the second a low bypass-ratio mixed-flow turbofan. 

4.3.1 Turboprop Engine 

Figure 21 illustrates the configuration selected for the 
sample turboprop test case. The engine has a single spool con- 
sisting of a two-stage centrifugal compressor and three-stage 
axial turbine. A reverse-flow burner and gearbox are also 
included . 

Figure 22 illustrates the WATE-S input necessary to estimate 
the weight and dimensions of this engine. The first: title card is 
followed by the configuration input contained in the $C namelist. 
This engine has 11 components and 11 thermodynamic flow stations, 
thus '‘NCOMP” and "NOSTAT" are set to 11. The component types are 
defined by the "JTYPE" array and the component connectivity Is 
specified by the "JCONF" array. 
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Figure 22. WATE-S Turboprop Sample Case Input. 


Following the second title oard^ the thecmodynamio input is 
specified in the $T namelist. The mass flow rate, total pressure 
and temperature, and fuel/air ratio are defined for each flow 
station. Note that the flow stations represent upstream component 
conditions. The fuel flow and brake shaft horsepower are entered 
in the "PERPP" array but this data is for informative purposes 
only and is not used by the WATE-S code. 

The aeromechanical input Is next, contained in the $W name- 
list. Default values are used for all but the "ILGNG", "INMEC", 
and "DESVAL” arrays. Basically, the default values set both the 
input and output units to English, turn on the weight estimation 
code and printer plot option, and provide a complete debug print- 
out. The length contributing vector "ILENG” does not contain DUCT 
6, PROP 10, and SHPT 11 since these components do not contribute 
to the total engine length. 

The "IWMEC values follow the "ILENG" array. The Inlet, com- 
ponent 1, is not entered since inlet calculations are not per- 
formed in WATE-S. Components 3 and 5 are specified as one-stage 
centrifugal compressors with a front frame associated with com- 
pressor 3. Component 4 is a crossover duct, IWMEC(2,4)«4 and com- 
ponent 6 is a burner. The HP turbine, component 7, derives its 
rpm from compressor 5, IWMEC(3,7) ■ 5y and is a three-stage axial 
turbine, IWMEC(5,7) = 3 and IWMEC(6,7) ■ 0. The propeller type is 
a Hamilton Standard double acting, advanced technology aluminum 
propeller, IWMEC{2,10) ■ 0. 

The "DESVAL" values follow the "IWMEC" inputs. "DEFAUL" val- 
ues are used for all but the compressors, burner, turbine, and 
propeller. Compressor 3 has an Inlet Mach number of 0.49, inlet 
hub/tip radius ratio of 0.356, and exit Mach number of 0.303. The 
number of blades is 17 and no splitters are used. DESVAL (4, 3) ■ 
17. The design point pressure-ratio tip-speed correlation Is 
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being modified 1.6 percent to match a known compreeeor rpm. 
DiJVVAT.(l3,.l) -- 1.0,16. liurner 6 has an airflow veloaity of 50. 
tl./s, rooidency time of 0.012 seoonda, and a mean diameter of 15.0 
Inch. 

Turbine 7 has an Inlet Mach number of 0.196, a moan work coef- 
ficient of 1.523, and an exit Mach number of 0,25, The tip solid- 
ity is 1.0 and the aspect ratios of the first and last stages are 
1.2 and 2,2, respectively. A constant hub radius design is speci- 
fied, DESVAL(8,7) - 3. The disk material Is a high strength 
auperalloy with a reference stress of 125,000 psl. 

The design cruise Mach number of propeller 10 is 0.,5 and the 
propeller has 4 blades. The activity factor is 130. and the tip 
speed is 600. ft/s. The propeller loading is 18. hp/ft^, 
DESVAL(5,10) - 18. 

Figure 23 illustrates the comp.'’ete debug output for the tur- 
boprop sample case. This includes the configuration, ^hermoc?/- 
namic, aeromechanical, and printer plot sections*, figures 2i ur.d 
25 illustrate the long and short form outputs of the Aero sohar- 
ical calculations of WATE-F. 
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Figure 23. WATE-S Turboprop Sample Case Debug Output. 
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Figure 23 . WATE-S Turboprop Sample Case Debug Output (Contd). 
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Figure 23. WATE-S Turboprop Sample Case Debug Output (Contd) 
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Figure 24. WATE-S Turboprop Sample Case Long Output, 
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Figure 2ii, WATE-S Turboprop Sample Case Short Output. 
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4.3,2 Turbofftn Engine 

Plaute 26 llXu.tt.tM the oondgutatlon ot the .ample tutbo- 
f.n teat oaae. The engine haa two apoola. The low-pteaaura apool 
conalata ot a two-atage axial tan, T-atage, ana two-atage axial 
low-pteauta turbine. The hlgh-pcee.ute .pool conalata of a 
alngle-ataga centrifugal oompteaaor arlven by a alngle-etage axial 
hi jh-preaeure turbine. The engine haa a revetae-flow burner and a 
ml xec “Compound exhaust system. 

Figure 27 llluatratea the Input requited to eatlmate the 
weight of thla engine. The tlret title card la followed by the 
configuration Input, namellat $C. Thla engine configuration la 
neatly Identical to the default and taqultea few ohangea. ■ ompo- 
nent 15 la a mixer Inetead ot a noaxle, and ahaft 16 oontulna 1 

fan, T-stage, and LP turbine. 

Following the aecond title card, the CberKOfynamlo Input La 
apeolfled In the $T namellat. The maaa flow rate, total pteeaure 
and temperature, and fuel/alr ratio ate apecltled for esch ata- 
tlon. The threat and fuel flow are entered In the "FERPF" array, 
but theae ate .ptlonal Inpute and are not requited tor execution, 
hirer, the mixer flow area., rATOOT,l,15, and PATOl.T(2,15, are 
requited and muat be Input for the mixer welgnt and dlmenalon cal- 

culations. 

The aeromeohanloal Inputa are next, and ate contained In the 
$W namellat. Default valuea have been uaed for moat of the 
“IMMEO” array since the eample tutbofan configuration la ao close 
to the default oonf Iguratlon. However, IMMEC(7,2) ’ > 

are set equal to 2 for the two-stage fan, and IMMEC(7,7) la se 
equal to 1 for the single-stage t-atage LP compressor. A two- 
stage LP turbine is Input, IMMEC(5,13) - 2 and component 15 la 

specified as a forced mixer. 
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Figure 27. WM’E-S Turbofan Sample Case Input. 



The "DESVAh” values follow the "IWMEC” inputs. In this case^ 
fttost '*DESVAL'‘ values are specified » the "DEPAUL" values used for 
1 . , -..w {.jomoonent 4), intercorapreasora duct (compo- 

nent 6 ) , turbine exit duct (component 14) , and shafts (compo- 
nents 16 and 17). The fan, (components 2 and 5) has an inlet Mach 
number of 0.56, inlet hub/tlp radius ratio of 0.4, and exit Mach 
number of 0.43. The tip solidity is 1,5 and the first and last 
stage aspect ratios are 2.0 and 1.5, respectively. The fan util- 
izes a constant tip radius design mode, DESVAL(12,2) «3.» and a 
3,5 percent scalar is applied to the design-point pressure-ratio 
corrected tip-speed correlation, DESVAL(13,2) = 0.965. The inter- 
compressor duct, component 8, has an inlet Mach number of 0.4 and 
a length-to-height ratio of 4.0. The inlet mean radius is deter- 
mined from the upstream component, LP compressor 7. The burner 
through-flow velocity is 55. ft/c and a burner residency time of 
0.008 seconds is specified. The inner diameter of the burner is 
specified as the outer diameter of component 9, the HP turbine 
(DESVAL(4,10) - -9.). This is how reverse-flow burners are 

handled when the mean diameter is not known. 

The HP turbine, component 11, has an inlet Mach number of 
0.16, turbine mean work coefficient of 1.24, and an exit Mach num- 
ber of 0.31. The tip solidity is 1.6 and the aspect ratio is 1.2. 
A disk reference stress of 125,000 psi is specified, DESVAL(7,11) 
= 125,000. The forced mixer has a length to hydraulic diameter 

ratio of 0.8 and has 8 lobes. 

Figure 28 illustrates the complete debug output for the tur- 
bofan sample case. This includes the configuration, thermody- 
namic, aeromechanical, and printer plot sections. 
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Figure 28. WATE-S Turbofan Sample Case Output (Contd) . 
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Figure 28. WATE-S Turbofan San^Jle Case Output (Contd) . 
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Figure 28. WATE-S Turbotan Sample Case Output (Contd) . 
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Figure 28. WATE-S Turbofan Sample Case Output (Contd) 


4 . 4 Program PlftgnostlQS 


The WATE-3 program contalne several error oheoks to aid the 
user in obtaining an acceptable engine design. A listing of these 
errors, and the subroutines the;' occur in, are included in 
Table XXV. None of these errors will cause the execution of 
WATE-S to stop. The subroutine the error occurred in may be ter- 
minated and subsequent component calculations can be in error. 

4 . 5 Program Structure 

The WATB-S code is written in FORTRAN IV and has been checked 
out on the IBM 370 computer. However, the code was developed at 
Garrett on a GDC 170 machine in an overlay structure. The code is 
written in single precision and requires no subroutines beyond 
those in the IBM FORTRAN IV manual. There is no character manipu- 
lation and only full word tests are used when testing BCD input. 

The main program is called WATES and in turn calls INPUT and 
WTBST as shown in Figure 29. Subroutine INPUT reads the configur- 
ation ($C) and thermodynamic ($T) data. Subroutine WTEST reads 
the aeromechanical ($W) data and calls the component routines. 
These component routines are independent of each other, and some 
use the same lower level routines as others. After all weights 
and dimensions have been estimated, ENGPLT is called to make the 
printer plot. 

The following variables in the common blocks may be refer- 
enced by a component weight estimating routine depending on the 
component type: DATOUT, WTP, TOPRES, TOTEMP, PAR, CORPLO, JCONP, 
JTYPE, NCOMP, NOSTAT. In no case is any value changed by the 
weight estimation code. 
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TABLt: XIV. WM’E-S ERHOtt MESSAGES 


JU^ROUTTNE FIC'JRB- 

«*1N PtOW f*TH OiONT 

noXPATH n*l OJtPNT CHEW lOSIC 


SU«Rf'UTI»*F niSK- 

*«**RIVER*«E DISK STRESS GREATER THAN ALLOMABII STRESS**R*R RRXXXXX. 


SURROUTINE OUCTl- 

' ouct "ts NOT converging" 


SUTRIUTINE CENm- 

STACF tl.PTSK STRESS EXCEEDS ALinVABlE. STRESS - XXXXXXX. AUCVAHE > XXXXXXX. 
STAGE TT GLADE STRESS EXCEEDS ALLOMAGLE* STRESS - XXXXXXX, ALLOWABLE > XXKXXXX. 
stage II blade HREOUENCY LESS THAN 2/REV 


SORRPUTINE CNECH- 

WARNING PILLOWING STAGE DESIGN LIMIT EXCEFDEO ***** 

BLADE ROOT STRESS 1? XXXXXXX. DES LIMIT IS XXXXXXX. 

AASTRESS IS TOO NIGH REDUCE SHAFT SPEED OR INCREASE EXIT NXCH NUMBER** 

**** WARNING FALLOWING STAGE DESIGN LIMIT EXCEEDED ***** 

STAGE HUATtP RATIO IS X.XX DES LIMIT IS X.XX 

**HUB TIP RATIO IS TOP HIGH REDUCE HUB TIP RATIO INPUT** 

**** warning FQUPWtNG STAGE DESIGN LIMIT EXCEEDED ***** 

ST*6E HUATIP RATIO IS X.XX DES LIMIT IS X.XX 

**HUB TIP RATIO IS TOO LOW INCREASE HUB TIP RATIO INPUT** 

**** WARNING following STAGE DESIGN LINIT EXCEEDED •**•« 

1ST stage PPATIO ts XX. XX DES LIMIT IS XX.XK 

**STAGE allowable pressure RATIO IS TOO HIGH REDUCE INPUT** 

**** WARNING FOLLOWING STAGE DESIGN UNIT EXCEEDED ***** 

LAST STAGE NACM NO IS X.XX DES LIMIT IS X.XX 

**LAST STAGE MACH NO IS TOO HIGH REDUCE STAGE EXIT MACH NO INPUT** 

**•* WARNING FOLLOWING STAGE DESIGN LIMIT EXCEEDED ***** 

STAGE BLADE HEIGHT IS XX.XX DES LIMIT IS XX.XX 

**STAGE BLADE HEIGHT IS TOO SMALL CHANGE OES ORR OR REDUCE H/T IN*UT** 

***« WARNING FOLLOWING STAGF OFSIGN LIMIT FXCEEOEO ***** 

S**GE WORN COEFFICIFNT IS XX.XX OES LIMIT IS XX.XX 

*P«TAGF WORK COEFFICIENT TO HIGH. ADO STAGES OR INCREASE H/T INPUT** 

COMPRESSOR II PRESSURE RATIO IS TOO HIGH 

COMPRESSOR It STAGE AND BLADE PARAMETERS* MEANINGLESS 
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TURBIME II STAGE AND SLACE PARAMETERS MEANINGLESS 



ORIGINAL tV.GS (£1 
OF POOR QlJAUTY 
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Baaed on the information in NCOMP and the JTVPE array, the 
proper component routine is called with the component number (I) 
as an argument. Each component is expected to fill WATE(X), 
AEiENG(I), TLENG(I), R0(1,I), R0(2,I), RI(1,I>, RI(2,1). Rotating 
components also fill RPMT{I). The shaft component fills DSHAPT(N) 


where N is the 

shaft count 

from the inside 

out. 

The array 

"CONVER” in 

common CONVER are conversion factors to 

convert English 

to SI units 

for output, or 

SI to English units for 

input. However, all internal calculations 
in English units. 

UNITS 

in WATE-S £.re performed 

LOCATION 

VALUE 

ENGLISH 

SI 

1 

2.5400 

inch 

cm 

2 

0.30480 

feet 

meter 

3 

0.45359 


Kg 

4 

0.092903 

ft^ 

2 

meter 

5 

0.027680 

lb /in^ 
m 

Kg/cm^ 

6 

0.68948 

Ibf/in^ 

2 

Newton/cm 

7 

4.8824 


Kg/m^ 

8 

0.55556 

*R 

®K 

9 

2326.0 

Btu/lbm 

Joule/Kg 

10 

0.074570 

hp 

KWatt 

11 

47.880 

Ibf/ft^ 

2 

Newton/m 


5,0 CONCLUSIONS AND RECOMMENDATIONS 

The WATE-S program can provide the level o£ aoouracy In pre- 
dicting weight and dimensions of both small and intermediate pro- 
pulsion gas turbine engines th-jt 1s required for preliminary 
engine cycle selection studies. On an absolute basis, the calcu- 
lated weights and dimensions have been demonstrated to be within 
±10 percent for six reference engines. On a relative basis, 
between engines of the same generic family and configuration, the 
accuracy is much better, estimated at less than ±2 percent. 

The flexibility of modeling the engine allows virtually any 
conventional or nonconvent I onal gas turbine engine to be analyzed 
given the thermodynamic performance of the engine components. 
WATE-S is viewed as an extension of the Boeing WATE-2 code, that 
is, the structure of the basic program was not altered. Hope- 
fully, this should allow users of the WATE-2 code to easily adapt 
WATE-S to their design system, if desired. 

The following Improvements to the program are recommended: 

(a) Exit Mach number should be used to size the blades for axial 
turbines instead of inlet Mach number as is currently in use. 
Exit Mach number is a more fundamental parameter in turbine 
design and the code should reflect this. 

(b) The heat exchanger calculation should be made to "close the 
loop" on fixed tube number. Currently, WATE-S estimates the 
pressure drops through the heat exchanger for the input num- 
ber of tubes. This could be compared against the thermody- 
namic data, and the number varied until the results agree. 
For nonpropulsion engines, e.g. auxiliary power units, 
stationary power units, and vehicular engines, the heat 
exchanger can be the largest and heaviest component in the 
engine, and hence, should be modeled as accurately as pos- 
sible. 
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(a) For ducts, burners, and duct burners where the ID wall is 
subjected to collapsing pressure rather than internal burst- 
ing pressure, the calculations outlined in Section 3.1.6 
should be included. Duct weights are typically a small frac- 
tion of total engine weight but a definite component accuracy 
Improvement could be obtained. In addition, the duct model 
could be Improved to radially match both the upstream and 
downstream components. 

(d) Blade mean solidity should be used to determine the blade 
count for axial components instead of the blade tip solidity. 
Blade mean solidity generally has more significance and con- 
sistency than does blade tip solidity. 


APPENDIX A 
I«I8T OP SYMBOLS 


2 

A area (om ) 

AF activity £aatoc 

AR aspect ratio 

C blade chord (om) 

C-D oonvergent>*divergent 

nozzle 

C/S solidity^ blade chord to 
spacing ratio 

D diameter (cm) 

F force (Newton) 

g gravitational constant 

Ki:M 
N-S"® 

GR gear ratio 

h height (cm) 

h/t hub/tip radius ratio 

H total enthalpy (Joole/kg) 

HP high-pressui e spool 

I rotational Inertia (cra^-kg) 

J 778.16 ft-lbf/Btu 

K factor for blade volume 

L length (cm) 

LP low-pressure spool 

Mn Mach number 

N number of elements 

P pressure (N/cm^) 

PTO power takeoff 

PW power (KWatt) 

R radius (cm) 

RPM revolutions per minute 

S blade spacing (cm) 

T temperature (®K) or 

torque (cm-N) 

t thickness (cm) 

TR blade taper ratio 


U wheel speed (m/s) 

1 

V volume (cm ) 

W weight (kg) 

T shear stress (N/cm^) 

\ turbine work 

coefficient 

P density (kg/cm^) 

u rotational velocity 

(rad/s) 

2 

<r normal stress (N/cm ) 

e heat exchanger 

effectiveness 

$ ratio of local total 

temperature to 
standard temperature 

d ratio of local total 

pressure to standard 
pressure 

iff compressor work 

coefficient 

Subscripts: 

h hub 

t tip 

B blade 

c case or corrected 

conditions 

S stator 

hw hardware 

Stg stage 

D disk 

SPL splitter 

2 engine inlet station 

m mean 

i inner 

o outer 
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